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In this paper, we study double charmonia production in Upsilon peaks, especially, a S-wave 
charmonium and a P-wave charmonium hd^PU), or a S-wave charmonium J/^ and the X(3940) 
and X(4160) within the nonrelativistic QCD (NRQCD) approach which is a powerful tool to realize 
the factorization of double charmonia production in electron-positron annihilation. The 
state T {nS) can provide an ideal laboratory for studying the properties of double-heavy quarkonium 
and also can separate the perturbative and nonperturbative parts due to the large heavy quark mass 
compared with the typical hadron scale Aqcd- Explanation of the X(3940) and X(4160) as the 
3^ So and 4^ So states, respectively, are compatible with the observed upper limits for the branching 
fractions of T(1S, 2S) —^ J/'if + X, where X = X(3940), X(4160) by the Belle Collaboration. The 
branching fractions of T(1S, 2S, 3S) —>■ ?7c + hc^Pi) are predicted to be around 10“®, which shall 
be tested in Belle-II experiments. 

PACS numbers 12.38.Bx, 13.25.Gv, 14.40.Pq 


I. INTRODUCTION 

The T(nS) below the threshold, as long-lived 

states with = 1 , have plenty of decay modes, 

the measurements which can be feasible in most current 
experimental techniques. Of particular interest to study 
those decay modes are the multigluon environment in the 
final products. Take the T meson for instance, its domi¬ 
nating decay mode is the decay into three gluons, which 
takes up 81.7% [1]. Another important mode is the ra¬ 
diative decay into two gluons, which possesses 2.2% [T]. 
Both of them lead to a multigluon intermediate state, and 
via short-distance hard interactions and long-distance 
soft effects, they furthermore transform into final measur¬ 
able color-singlet states including conventional hadrons 
predicted by the constituent quark model, and exotic 
states beyond the constituent quark model such as hy¬ 
brids, multiquark states, hadron molecules, and Glueball. 
These = 1 states with highly narrow decay width, 
also including charmonia ijj{nS) below the thresh¬ 

old, have natural advantages to hunt for exotic states 
since there are a great many events for these vector states 
that strongly couple to multigluon. Lattice QCD also in¬ 
dicates there are a large possibility for Glueball product 
in vector heavy quarkonium decays [2 [3]. To more pre¬ 
cisely extract the signals of the exotic states from the con¬ 
ventional spectrum, however, one need to study clearly 
the conventional hadrons production in those multigluon 
product processes. 

The large bottom quark mass provides a hard scale 
to separate the contributions of the strong interactions 
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into short-distance and long-distance effects. This kind 
of separation is wellperformed in the nonrelativistic QCD 
(NRQCD) g, where the amplitude is expressed by a sum 
of products of non-perturbative NRQCD long-distance 
matrix elements (LDMEs) and the corresponding short- 
distance Wilson coefficients. A proof of the NRQCD fac¬ 
torization in the exclusive double-charmonium produc¬ 
tion in an e"*" e~ annihilation, and a quarkonium and a 
light meson production in B decays was recently given 
in Refs. In this paper, we employ NRQCD factor¬ 

ization and study the double-charmonium production in 
the T{nS) decays. 

Double-charmonium production has attracted a large 
amount of attention among theorists and experimental¬ 
ists in the last ten years, especially the large discrep¬ 
ancy between theory [THin] and data [ni [la in the 
B factory, that is later solved through QCD correc¬ 
tions HMD- In recent years, some exclusive decay chan¬ 
nels of bottomonium into double charmonia have been 
investigated intensively, e.g., T —>• J/if + r]c [TM2fi] . 

r]b -t J/ii + J/if [21]^, XhJ -t J/'4^ + J/i’ [25H28], 

T — >■ J/ lit + XcJ US]- The branching fractions of these 
channels are predicted to be order from 10“® to 10“® in 
the above literatures. In experimental aspects, the Belle 
Collaboration has recently measured the branching frac¬ 
tions for the T(1S', 2S) decay into double charmonia [50] . 
Except the known charmonia, the Belle Collaboration 
has also observed the A(3940) |3T] and A(4160) [32] sig¬ 
nals through the spectrum of mass recoiling against the 
J/ijj in the process e+ -I- e“ —>■ J/V' + X. Up to now, to 
interpret the constructions of the X(3940) and A(4160) 
states is still a challenging issue due to lack of enough 
data. In this paper, we take the A(3940) as the 3^S'o cc 
state and the A(4160) as the 4^5'o cc states, and attempt 
to find out whether or not to explain the results observed 
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in the Belle experiment. We will give the amplitudes of 
both T{nS) —>■ hd^Pi) + rjc and T{nS) —)• J/'0 + X, 
where X = X(3940), X(4160). 

The paper is organized as the following. In Sec. |T^ 
we present the NRQCD factorization calculation formu¬ 
lae. The amplitudes of the channels in question are given 
accordingly. In Sec. |III[ the branching fractions are pre¬ 
dicted and the comparisons with data are also presented. 
We summarize and conclude in the end section. 


II. FACTORIZATION FORMULAE 


V between the heavy quark and heavy antiquark with 
the mass mq in the heavy quarkonium H. Im/„(/rA) is 
the corresponding short-distance coefficient, which can 
be calculated order by order in perturbative theory. 

According to the order counting rules, the lowest-order 
NRQCD operators which contributes to the above pro¬ 
cesses we concerned are [3] 

(3) 

C)(3S'W) = il^'fcrx- (4) 

OCPW) = (5) 


A. NRQCD approach 


The NRQCD Lagrangian includes three parts, i.e. the 
heavy piece which describes a Schrddinger held for each 
heavy quark or antiquark, the light piece which de¬ 
scribes ordinary QCD without heavy quarks and heavy 
antiquarks, and the correction piece which ensures that 
NRQCD reproduces QCD results. It can be written as 
follows [3] 

if 

>Cnrqcd = ^ ^ ^ 

+ • (D X 5,E - gsE X D)^ 

-I- (^/> —>■ zCT^x*, —>• — A^) -I- £iight, (1) 

where ip and x denote the Pauli spinor held that an¬ 
nihilates a heavy quark and creates a heavy antiquark, 
respectively. Dt and D are the time and space compo¬ 
nent of the gauge-covariant derivative . A® = 
and S® = ied'cc'jfc are the electric and magnetic color 
components of the gluon held strength tensor G^®®. The 
replacement in the third line implies that one can ob¬ 
tain the corresponding heavy antiquark bilinear sectors 
via the charge conjugation transformation. Gught rep¬ 
resents the Lagrangian for the light quarks and gluons. 
The coefficients ce, cf, and cq have perturbative series 
in powers of the strong coupling Us , which can be written 
as Ci = 1 -|- 0{as). 

Before we write the NRQCD factorization formulae 
for the amplitudes of both T(nS') — hc{^Pi) + gc and 
T{nS) —>■ J/ip + X, where X = X(3940), X(4160), we 
introduce the corresponding NRQCD LDMEs. Factor¬ 
ization of the inclusive annihilation decay width of heavy 
quarkonium can be written as 

m) =Y,‘^h^{H\Or.{^JiE)\H) , (2) 

TTlr^ 

71 Q 

where LDMEs {H\On{g-A)\H) involve nonperturbative 
effects and are well-organized by the relative velocity 


The corresponding matrix elements of the operators 
sandwiched by meson states are 

= (iL|G(2«+iLW)|iL). (6) 

B. The amplitude for T{nS) + rjc 


In this subsection, we calculate the amplitude for 
T{nS) hc{^Pi) + rjc- Before performing the calcu¬ 
lation, we introduce an equivalent method mUH], i.e. 
the covariant projection method, rather than the direct 
matching method. 

The Dirac spinors for the heavy quark with momentum 
Pi and heavy antiquark with momentum p2 in quarko¬ 
nium can be written as 


uq(gi, A) 


VQiP2,>^) 


El -I- mg 


2Ei 


E2 + mg 


2 E 2 


Ei+niQ J 


3-V2 


E'2+mQ 



(7) 

( 8 ) 


where Ei are the corresponding energy of heavy quark 
and heavy antiquark, which satisfy Ei = E 2 = E. We 
introduce q as half relative momentum between the heavy 
quark and heavy antiquark with pn • q = 0, where pn = 
Pi + P 2 - In the rest frame of quarkonium, we have E = 

ijrriQ — q^. is the two-component Pauli spinor and A 
is the polarization parameter. Using the above formula 
for Dirac spinors, it is straightforward to get the covariant 
form for the spin-singlet and spin-triplet combinations of 
spinor bilinearities. The projections are 


ns(<7) = 


UQipi,Xl)vQip2,X2){^Xl^X2\SS^) (g) 

+ 2i? 


Ai ,X2 


+ ^- mg) (g , 


(9) 
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where = 7^ for the spin-singlet combination, while 
Fg^i = ^{ph) = 'y^Sfj.iPH) for the spin-triplet com¬ 
bination with the polarization vector £^(jph)^ and the 
spin-singlet projection no(g) and the spin-triplet projec¬ 
tion ni(q) are defined accordingly. Ic is the unit ma¬ 
trix in the fundamental representation of the color SU(3) 
group. Besides, one needs to note that the state |iF(p)) 
in NRQCD has the standard nonrelativistic normaliza¬ 
tion: {H{p')\H{p)) = (27r)^(5^(p — p'), while an addi¬ 
tional factor 2Ep is included in the state normalization 
in QCD where {H{p')\H{p)) = 2Ep{2Tr)^6^{p — p'). 

To obtain amplitudes of P-wave state production, one 
can do the Taylor expansion of the amplitudes in powers 
of 


A{q) 


^(0)-h 


dA{q) 

dqt^ 


q=o r 


1 d^A{q) 
”^2! dqf^dq’' 


3=0 +- 


( 10 ) 




FIG. 1: Typical Feynman diagrams that contribute to the 
process T{nS) —>■ hcCPi) -I- rjc. 


The amplitude for T ^ he + rjc can be expressed as 
two independent terms 


MiT^he + Pc) = El e{pr) ■ £*{phj 

+FoPhe ■ £{pr)pr ■ AiPhJ- 

( 11 ) 


In the nonrelativistic QCD approach, the amplitude of 
T —>■ he + 7 c at leading-order in v can be factorized as 


where 


and 


M{T^hc + Pc) = '^N{lL^,00\l,J^){0\x’cr'il;\T){hc\'tpH-^^)x\0){Vc\'fk'’x\0) 

Lz 

X Tr [71'^" (0)no (0) + 71“^ (O)n^ {pn ^, L,)e" {pr ), 


i (N^ - 4) (N^ - 1) 7“.7’^.(2zTOfc -f - I^i + 

128 v^ 7 rz 2 m 37 Vc'^^ -I- fc2)^ ql (hi - hi - {{qi - kif - z'^ml) 


■Tr[(2mb -I- ^r)-Y-7^-(mb + i 

1 - 

- ^).7“.(mb + - |(5i - % 

+ ^)y] 

( 

(gi - -ml^ 

( 

{qi-ki-'F^ -ml^ 



^ Tr[(2m{, -h;iT)-7''-7“-(wb - ^ + ^)-7^-(wb - ji + ^)-7^] 

((gi - mi) ((hi - mij 

^ TT[{2mb+^T)-YA^-imb - ^2 - - %^ + ^).7'^.(m{, - ^2 - ^ + ^)-7°] ' 

((gi + k2 + ^-^i- mi) ((fc2 


n^(o) 


5no(g) I 
dqf^ 




( 12 ) 


(13) 


(14) 


N = y/2m-Yy/2mh^yj2'mp^l{ 2 Nc)^A is from NRQCD 
operators and states normalization factor, and the 
vacuum-saturation approximation for LDMEs is also 
used: {H\0^\H) - (iF|V^t^;^|o)(0|xt^„,/,|i7) ^ith 

For convenience, we have de¬ 
fined z = rric/mb, and the momentum ki of charm 
quark(antiquark) in the hadron can be written as 


hi = ^ -I- g and k 2 = ^ — q. 

For the charmonium state the summation over the 
quark spins and orbital momentum projections results in 

^(lL„00|l,JQe*'^(p;,^,L,) = £*^(p,„,J,). (15) 

L:, 
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The summation over the polarization of he is 

1 f-^ ly 

J^=-l 

Next let us give the explicit amplitude of T —>■ he + 
Tje- For pure hadronic decay, typical Feynman dia¬ 
grams are depicted in Fig. There are twelve dia¬ 
grams contributing to the process, and another nine di¬ 
agrams can be obtained by reversing direction of each 
quark line one by one. According to the power-counting 
rules [1ES], the matrix element scales 

as :y3+2i+2E+4M^ where S and L are the spin and or¬ 
bital angular momentum quantum number for the QQ 


pair, and E and M are the minimum number of chromo¬ 
electric and chromo-magnetic transitions for the QQ 
pair from the dominant Fock state of H to the state 
QQ (^^~^^Consider the higher-order in v, the non¬ 
trivial next-to-leading order contribution for the above 


process are |56 (^S'P))t 


\bb{^PyUi,2))r 


cc(^S'o®'))h<, + \ccysy))r,,, and 




\bbyPjlo,i, 2 ))r |cc(ip|^'))/,^-f|cc(ipf')),,^, where the 
first one has a relative suppression uF, while the last two 
have a relative suppression compared to the color- 

singlet contribution. 

The diagrams in Fig. [^together with other nine topo¬ 
logical diagrams give the corresponding coefficients 


F’n = -- 




32:3(1 — 4z^)^(3z2 — 


x{l920z^° - 4256z® -f 432z® -k 324^^ - 372^ + 3 -h (-1792z® -f 640z® - 20z^ 
-13z^ + 3)6i -f 12 z2(80z® - 44z® - 22z* + 132^ - 1)62 + 22^(322 - l)(152z4 
-f38z2 - 1)63 -h 4 z2(340z® - 1472-^ -f 30^^ - 7)&4 + {5z^{-520z^ + 16z® 
-f96z^ - 58^2 -f 17) - 3)65 + 6z2(3z 2 _ i)(i6z® -h 202"^ -f lOz^ - l)ci 

- 12 ^ 2 ( 3^2 _ g ^2 


(17) 




2ttCf i2NeCF - 3) gf 

3^3 ( 42;2 _ 1 ) 


(01 yl'erV' IT) (he I( - ^ ^) X 10) ( 77 c I X 10) 


X { 112^4 -k 82^ -f 3 -f ( 16 z^ -k 8 z 2 -f 3 ) 5 i -k 12 z'^{ 4 z^ - 1)62 
- 2 z ^{ 20 z '^ + 1)63 -f 242^(2^ _ 1)64 + (-482^ -f 302 ^ - 3)65 

- 62^(822 - l)ci - 1222(52^ - 2 )c 2 }, 


(18) 


where the one-loop master integrals bi and Ci are pre¬ 
sented in the appendix. 


C. The amplitude for T{nS) —^ J/ip + X{M) 

Considering the Belle Collaboration has observed the 
channels of T(n5') —>■ J/ip + X where X can be ei¬ 
ther one of the A (3940) and A(4160), we go to study 
these channels. The A(3940) [3T] and A(4160) [32] are 
first observed by the Belle Collaboration in the spec¬ 
trum of mass recoiling against the J/if) in the process 
e"*" -I- e“ —>■ J/t/j -I- X. Up to now, we have no more 
internal information except the fact that both of them 
shall have a charm quark-antiquark pair because of their 
hadronic decays into DD or D* D. There are some dif¬ 
ferent schemes to interpolate the nature of the A(3940) 
and A(4160), e.g. a hybrid charmonium |34j . a molecular 


charmonium [35j , and a pure charmonium state |361l39j . 

In this paper, we take in a pure charmonium expla¬ 
nation for these states, and adopt the 33 50 and 4^50 
states to interpolate the A(3940) and A(4160), respec¬ 
tively. Then we just calculate the process of T(n5) —>■ 
J/ijj + cc(335o, 4350 ). The pioneer works of T(n5) —>■ 
J/iIj + T]c have been performed in Ref. (TS] [T9j. We have 
also calculated the amplitudes, which are in agreement 
with Ref. [15] . 

Following the same procedure, the amplitude of T —>• 
J/t/j-I-A( 3940) at leading-order in v can be factorized as 

M(T^ J/V'-f A(3940)) 

= N'{0\x^(Tip\T)){J/'ip\'ip'f(Tx\0) 

x{X{3M0)\yhmMy0)U^{0)j/^] 

'y£*^{pj/i,)e''{.Pr), (19) 

where N' = •\/3-\/2mxA/2mj/.0-\/2mjs:/(27Vc)3/2 and 
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^‘^(0) are identical to Eq. at q = 0. 

The amplitude for T —>■ J/^/> + X(3940) can be written 
as 

J/V' + X(3940)) 

= iFAef,^ap£^{PT)e*''{pj/^)prPh^. ( 20 ) 


The diagrams contributing to T —)• J/^/; + X(3940) are 
analogous to T —>■ he + rye, which can be got by replacing 
he to J/ij] and to X(3940) in Fig. We obtain the 
coefficient 


z‘^ {4:Z^ — 1) 

x{-16z^ + 12z^ - 1 + (4^2 _ l)5i - 4 z2 {2z^ - l) &2 - 

X - Az^ + (Sz^ - 2z^ + 1)65 + 2z‘^ {2z‘^ - 1) Cl + Sz^Cs}. (21) 


Because we treat the X(3940) and X(4160) as the 3^S'o 
and 4^S'o states separately, the coefficient for T —>• 
J/ip + X(4160) can also be obtained by just replacing 
the LOME (0|xti^|x(3940)) to (0|xti/,|x(4160)). 


III. PHENOMENOLOGICAL DISCUSSIONS 


Experimental data by the Belle Collaboration indicate 
that the cross section of e+ + e“ —)■ J/^tp + X where 
X = A(3940), A(4160) is not trivial [311 [32]. One has 
to resort to certain schemes to describe the nature of the 
A(3940) and A(4160), and different treatments deter¬ 
mine distinctive results. A pure charmonium explana¬ 
tion for both the A(3940) and A(4160) are popular in 
current literatures |37H39] . Of course, the measurement 
with more data will tell us the nature inside the A(3940) 
and A(4160). Except for the A(3940) and A(4160), mea¬ 
surements of double charmonia production with a S-wave 
charmonium and a P-wave charmonium are also feasible 
in current Belle experiment [30] . 


In the above section, we have given the NRQCD factor¬ 
ization formulae and calculated the corresponding short- 
distance coefficients. Next we input the value of LDMEs 
and other fundamental parameters. In this paper, the 
A (3940) is treated as the pd^S), and the A(4160) is 
treated as the pc{4:S). For the LDMEs of pc{nS), we can 
get relation the relationship of y)’(nS') using the heavy 
quark spin symmetry. In the leading-order of heavy 
quark relative velocity v, they are identical. So the 
LDMEs for the A(3940) and A(4160) can be obtained 
correspondingly. For the LDMEs of ip{nS), we can obtain 
the value from the measured electric widths [4| [401 [4T| , 
i.e. 


(0(35™))^/^ 


T{J/d ^ e+e-) 

87ra2e2 (1 — 4Q;sC'F/7r) 


( 22 ) 


Using the above expression, the corresponding values of 


both LO and NLO for the LDMEs are given in Tab. |l| 
The parameters we adopted are [T] 

my = 9.4603GeV, TOy( 2 S) = 10.023GeV, 

”^r( 3 S) = 10.355GeV, my( 45 ) = 10.579GeV, 
TOy( 5 S) = 10.876GeV, = 2.984GeV, 

^r;c( 2 S) = 3.639GeV, mj/^ = 3.097GeV, 
ml( 2 S) = 3.686GeV, ruh, = 3.525GeV, 
w^(404o) = 4.039GeV, m^( 44 i 5 ) = 4.421GeV, 
r(T) = 54.02keV, r(T(2S')) = 31.98keV, 
r(T(3S')) = 20.32keV, ree(T) = 1.285keV, 
ree(T(2S')) = 0.612keV, ree(T(3S')) = 0.443keV, 
ree(T(45')) = 0.272keV, ree(T(5S')) = 0.31keV, 
TeeiJ/d) = 5.55keV, ree(V(25')) = 2.36keV, 

Tee(^/^(4040)) = 0.86keV, ree(V'(4415)) = 0.58keV. 

For the mass of the A(3940) and A(4160), we have 
wx( 394 o) = 3.942GeV and mx( 4 i 60 ) = 4.156GeV. The 
heavy quark mass is adopted as rric = 1.5GeV and 
mb = 4.8GeV [HI |33]. The strong coupling constant 
is set at the Z-boson point with asijnz) = 0.1185 where 
mz = 91.1876GeV [I], so we can run the coupling to 
other points. For example, Q;s(4.8GeV) = 0.2178 at two- 
loop evolution with active flavor n/ = 5. 

For the LDMEs of Pc{nS), we can obtain the following 
approximation using the heavy quark spin symmetry 

(0(1.5W)),,(„S) (0('^i'’))^(n5). (23) 

To take V'(4040) and ?/!(4415) as '0(45') and 0(55) respec¬ 
tively [HI, we have 

(0(l5W))^,(35) = (O('5 W))v-(4040) 

= 0.1083;°;°3“(GgV)3, 

(0(l5W))^^(45) = (0('5 W))v-(4415) 

= O.O875;°;°2?4(GeV)0 


(24) 
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TABLE I: The value of LDMEs for i>{nS) and T(nS') ex¬ 
tracted from the electric widths at LO and NLO, in compari¬ 
son with the results from B-T potential model by Buchmiiller 
and Tye |45[ I46| . The strong coupling constant is evaluated 
at the scale 2mQ, and varying the scale from mg (2 to the 
meson mass we can get the corresponding uncertainties. 


(GeV)^ LO NLO B-T model 
0.2341 0.4107t°;“78 0.3867 

(C>(3S'f'))^(2S) 0.1411 0.2475t°°®®i 0.2526 

(C>(3S'f'))^(3S) 0.0617 0.1083t^;;;?®^ 0.2172 

(0(35W))^(4S) 0.0499 0.0875t°:°??^ 

2.025 2.925155;^® 3.093 

(0(®S'f’))T(2S) 1.082 1.563t°;J®:^ 1.544 

(0(®5W))t( 3S) 0.8361 1.208t°:“® 1-181 

(C>(®S'f'))T(45) 0.6340 0.9158t2jj^°7 1.025 

{0(®gM))Y(5S) 0.6413 0.9262tg °fJ 0.9339 


For the LDMEs of we adopt the result of B-T po¬ 
tential model in Refs. 05] 

{OCPl^^))h, = 0.1074(GeV)®. (25) 

The decay width of T —> he + rye cs-n be written as: 

r(T ^ he + rye) = -^|M(T^he + r?e)P, (26) 

oTTTTZ'Y' 

where 

I I _ y/ (to| - {iTlrj^ - ITlhJ 2) (to| - {mh, + TO^J 2) 

'P' " 2 mr 

is the momentum modulus of final charmonia in the T 
meson rest frame, and Ai = is the amplitude 

with average polarization. 

In the end, we give the branching fractions of the 
processes we concerned, which are presented in Tab. [l^ 
The Belle data are also listed for comparison. The re¬ 
sults we presented can explain the experimental upper 
limits for the processes in question. And we also pre¬ 
dict the branching fractions of T(rrS') —>■ he + rje and 
T(3S') J/Ip + X{M), which have the potential to be 
observed in current experiment. Note that in Ref. m, 
Jia has investigated the process T —>■ J/ip + ijc includ¬ 
ing both the three-gluon mode and the electromagnetic 
radiative mode contributions. However, the electromag¬ 
netic radiative mode contributions are small as the au¬ 
thor has pointed out. 


IV. CONCLUSION 


In this paper, we have calculated the branching frac¬ 
tions of T{nS) he + rjc and T(nS') —)■ J/ip + X with 


TABLE II: Branching fractions (10“®) for T{nS) hadronic 
decays into double charmonia, where the first column uncer¬ 
tainty comes from the scale running from 2mc to 2mi, and the 
second column uncertainty is from the heavy quark mass with 
rric = 1.5 ± O.lGeV and m;, = 4.8 ± O.lGeV. For comparison, 
the results Bjia, are from Ref. m, and the data are from the 
Belle experiment m- 


Channels 
T J/p) + rjc 
r ^ j/pj + nc{2S) 
r ^ j/ip + A (3940) 

T J/V> -t A(4160) 

T —> Tjc he 
T ^ xP{2S) + rje 
r ^ J/iP + A(4160) 
T(25*) —y Tjc he 
T(35') J/lP + Tie 
r{3S)^pji2S) + Vc{2S) 
T ip{2S) + X(3940) 
T ^ V’(2S') -b A(4160) 
T(25') J/r/> -b rie 

T(25) ^ J/V> + Vci2S) 
T{2S) J/ip + X{3M0) 
T(2S) ^ J/'iP + v42S) 
T(3S') J/ip + A (3940) 

T(3S') ^ J/^A-b A(4160) 

T(35*) —y Tje he 


I^NRQCD_Hjia I^Exp 


Q Q9+5.77+0.59 
'^•^^-2.62-0.53 

3.9+2®3 

<2.2 

1 'y<y+2.61+0.27 
' -1.18-0.24 

2.0-pii 

<2.2 

A f-o+1.00+0.10 
U.UO_Q 45_Q Qg 

- 

<5.4 

0 4Q+0.74+0.07 
^•^^-0.33-0.07 

- 

<5.4 

■t qo+ 1-96+0.34 
i.OO_o.89_o.26 

- 

- 

Q -1 g+3.21+0.33 
^•-10-146-0.29 


<3.6 

r, 40+0.64+0.06 
'-'•^'>-0.29-0.06 

- 

<2.0 

. Tq+ 1.65+0.29 
i-i^_0.75-0.22 

- 

- 

Q or.+5.59+0.57 
o-OU_2.54_o.5i 

- 

- 

Q7+1.42+0.14 

-0.64-0.13 


<3.2 

^ q<:‘+0.54+0.05 
U.OO_o.24-0.05 

- 

<2.9 

A q^+ 0.39+0.04 
U.ZO_oj^g_0.o4 

- 

<2.9 

O 97+4.82+0.49 
' -2.19-0.44 


<5.4 

1 oQ+2.05+0.21 
-L-oU_o.93_o.i9 

1 0+2.1 
-‘-•'^-0.9 

<2.5 

^ re+0.86+0.09 
U.OO_o 39 _o.o8 

- 

<2.0 

n ^4+2.41+0.25 
-‘-•^^-1.09-0.22 

- 

- 

A f-D + 1.00 + 0.10 

o-L10_o.46_o.09 

- 

- 

A r 1+0.75+0.08 
Ll-Oi-o.34-o.07 

- 

- 

1 oi +1.92+0.33 
i-oi-o.87-o.25 

- 

- 


A = A(3940) or A = A(4160) using the NRQCD fac¬ 
torization approach. For the Upsilon peaks below the 
threshold, the branching fractions of the processes 
in question are around I0“® — 10~^; thus, these channels 
have the potential to be measured in the Belle-II exper¬ 
iment, where the events of T, T(2S), and T(3S') will go 
to 1.8 X 10^^, 7.0 X 10^°, and 3.7 x 10^°, respectively 07]. 
At Upsilon peaks higher than the threshold, the 

branching fractions will be suppressed by a factor around 
10^, because the decay width of the corresponding T(nS') 
meson increases greatly up to dozens of MeV. The Up¬ 
silon family, also existing in the J/pj family, provides a 
natural laboratory for multigluon products and indicates 
the potential existing space for the exotic states beyond 
the constituent quark model. NRQCD provides an ef¬ 
fective factorization formula for the double-charmonium 
production, where the short-distance Wilson coefficients 
and the LDMEs are well-factorized, since the bottom 
quark mass is a large scale compared with the hadron 
scale Aqcd- But we should also notice that the factor¬ 
ization to all orders of the strong coupling as and all or¬ 
ders of the quark relative velocity v in Upsilon hadronic 
decays into double charmonia is still needing further in¬ 
vestigation. 
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To understand the data carefully, one should clarify 
the nature of distinct hadron states. Up until now, 
the physics information of the X(3940) and X(4160) 
is lacking except for their masses and poorly studied 
decay modes. The explanations of the X(3940) and 
X(4160) as the 3^S'o and charmonium states, re¬ 
spectively can explain the experimental upper limit for 
T{nS) —)■ J/ip + X{M). The branching fractions are 
around 10“^ when we treat the X(3940) and X(4160) as 
a pure charmonium, and future measurement with more 
data by Belle Collaboration will be able to clarify it is or 
not. 
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Appendix 


In our calculation, the package FeynCalc [48] is used 
to generate the Feynman amplitudes, LoopTools [35] is 
used to calculate the integrals, FIRE m and Apart m 
are employed to reduce Feynman integrals to Master in¬ 
tegrals. The scalar Passarino-Veltman Master integrals 
Bi and Ci are defined in Refs. |3S1IS11IS3|, and bi and a 
are related to them with the following identity: bi = Bi, 


Cilml. 



Bi = 

Bo 


B2 = 

Bo 


Bs = 

Bo 

(mg, 0,0) , 

Bi = 

Bo 

{ml,ml,ml) , 

B5 = 

Bo 

[z m^) 

Cl = 

Co 

( 2 2 2 2 2 

C2 = 

Co 

( 2 2 2 2 2 
\m^,z m^,z 
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